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Summary. We introduce the seminal developments in the theory and experiments 
' of electromagnetic probes for the study of the dynamics of relativistic heavy ion 

^Sj I collisions and quark gluon plasma. 



1 Introduction 

O 

^ ' Collision of heavy nuclei at relativistic energies is expected to lead to formation 

I of a deconfined state of matter known as Quark-Gluon Plasma (QGP) [1], 

where quarks and gluons are the effective degrees of freedom rather than 
nucleons or hadrons [2, 3]. It is now well accepted that a few microseconds 
^ ' after the 'Big Bang' the whole universe was in the state of QGP [4]. 

, Several experiments performed at the Super Proton Synchrotron (SPS) at 

CERN and Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory, New York, have provided a significant evidence of the formation of 
this novel state of matter. A giant accelerator known as Large Hadron Collider 
(LHC) at CERN will be in operation very soon and will provide many new 
insights about the properties of QGP and the theory of strong interactions. 
\ Heavy ion collisions at relativistic energies produce extremely high temper- 

atures and energy densities within a very small volume. As a result, quarks 
and gluons (also known as partons) no longer remain confined within the 
' nucleonic volume and create a deconfined state of partons due to multiple 

scatterings and production of secondaries due to gluon multiplications. The 
system (or the fireball) may reach a state of local thermal equilibrium. It cools 
by expansion and below a certain critical temperature (Tc ~ 180 MeV) or en- 
ergy density, the partons are confined to form hadrons and the system reaches 
a hadronic state. It may undergo a further expansion and cooling before the 
freeze-out takes place. 

Radiation of photons and dilcptons has been proposed as the most promis- 
ing and efficient tool to characterize the initial state of heavy ion collisions. 
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Unlike hadrons, which are emitted from the freeze-out surface after undergo- 
ing intense re-scatterings, photons come out from each and every phase of the 
expanding fireball. Being electromagnetic in nature, they interact only weakly 
and their mean free path is larger than the typical system size {'^ 10 fm). As a 
result once produced, they do not suffer further interaction with the medium 
(a <C as) and carry undistorted information about the circumstances of their 
production to the detector [5]. 

Initially, photons (real as well as virtual) , were studied in order to get only 
the temperature of the plasma. Several other possibilities, e.g., (i) evolution 
of the system size by intensity interferometry [6] , (ii) momentum anisotropy 
of the initial partons [7] as well as formation time of quark-gluon plasma [8] 
using elliptic flow of thermal photons, (iii) an accurate check on jet quenching 
and other aspects of the collision dynamics by photons due to passage of high 
energy jets through plasma [9] etc. have come to the fore. Of-course dileptons 
are considered as the most reliable messengers of the medium modification of 
vector mesons [10]. 
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Fig. 1. Schematic diagram of different sources of pliotons and their relative pr 
spectra. 



2 Sources of photons 

In order to proceed, it is useful to identify various sources of photons from 
relativistic heavy ion collisions. Their production is a result of convolution 
of the emissions from the entire history of the nuclear collision. Photons are 
emitted from the pre-equilibrium stage, from QGP phase, from hadronic phase 
and also from the decay of hadrons produced at the time of freeze-out. An 
ideal situation would ensue if the contributions from different stages domi- 
nate different parts of the pr spectrum. A schematic diagram of the different 
sources of photons and their slopes is shown in Fig. 1. We, thus need rates 
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and models to study the evolution from different sources. Hydrodynamics, 
Cascade, Fire-ball, Cascade+Hydrodynamics are the vastly used models for 
this purpose. 

2.1 Direct photons 

The term 'direct photons' stands for the photons which emerge directly from 
a particle collision. In a heavy ion collision experiment, the detector captures 
all the emitted photons including those from decay of final state hadrons. The 
resultant spectrum is the inclusive photon spectrum. However, more than 90% 
of the photons in this spectrum are from hadron decay. One can subdivide this 
broad category of 'direct photons' into 'prompt', pre-equilibrium, 'thermal' 
(from QGP as well as hadronic phase) and 'jet conversion' depending on 
their origin. Before we come to the different sub-categories of direct photons, 
we start our discussion with decay photons and their subtraction from the 
inclusive photon spectrum. 

2.2 Decay Photons 

As mentioned earlier, most of the decay photons are from 2-7 decay of 7r° and 

rj mesons, lo, rj' etc. also contribute to the decay photon spectrum, marginally. 
Subtraction of the decay background from inclusive photon spectrum is a very 
challenging task. WA98 Collaboration [11] used the subtraction method using 
invariant mass analysis for decay background and later PHENIX Collabora- 
tion [12] has developed this method to a much higher level of sophistication. 

Invciricint mass cinalysis 

Extraction of direct photon spectrum from the inclusive photon spectrum is 

done using invariant mass analysis as the primary step. First of all, all the 
detected photons are listed on an event by event basis. Then by selecting two 
photons randomly from an event, invariant mass of the pair is calculated. If 
El, E2 are the energies and pi, p2 arc the three momentum of the photons 
respectively, then the invariant mass of the pair is, 

M^^=[(Ei+E2)2-(pi+p2)2]l/2. (1) 

If the value of M-y^ is close to m^o , it is assumed that they arc the decay 
products from the same pion and a 7r° spectrum is obtained. Similarly the rj 
spectrum is obtained. These two spectra are then used to determine the decay 
photon spectrum using kinematics and the subtraction of this spectrum from 
the inclusive spectrum gives the direct photon spectrum. 

However, in an event having N number of photons, the total number of 
photon pairs that can be formed is ^C-z- For large values of N, this number 
gets very large. Thus, there is a very high probability of getting a pair of 
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Fig. 2. Invariant mass distributions of pairs of electromagnetic clusters passing 
photon selection cuts for pair transverse momenta satisfying 3.0 < 3.5 GeV. Top 
panels: m^-, distributions in Au + Au events compared to a normalized mixed- 
event sample representing the combinatoric background. Bottom panels: The m-,^ 
distributions after subtraction of the combinatoric background for 10-20% and 60- 
70% centrality bins [13]. 



photons having an invariant mass m^n-o or m^, which are not decay product of 
the same pion or eta meson. 

Several other issues also need to be accounted for: (i) the detector resolu- 
tion is finite, (ii) the opening angle between the decay photons can get very 
small, specially if the energy of the pion is large, (iii) the photons may not 
deposit all their energy in the detector, and (iv) one of the photons may be 
outside the coverage of the detector. 

Mixed event analysis 

Thus, the major problem of invariant mass analysis is that, the accidental 
(false) photon pairs can also give rise to pion mass and it is not possible to 
distinguish them from the correlated pairs. To overcome this problem, a mixed 
event analysis [13] procedure has been used successfully. The basic idea of 
mixed event technique is to compare particle spectrum from one event to the 
result for particle combinations from different events, which are a priori not 
correlated. As a first step, properly normalized mixed events are constructed 
by randomly sampling photons from different events. The difference of the 
invariant mass spectra (see Fig. 2) of the real events and the mixed events then 
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gives the pion and rj distributions. Once again, the decay photon spectrum is 
subtracted from the inclusive photon spectrum to get the direct photons. 



Internal conversion and tagging of decay photons 

An alternative approach of separating direct photons from decay background 
is by measuring the 'quasi-real' virtual photons which appear as low mass 
electron positron pair. It is assumed that, any source of real photons also 
produces low mass virtual photons which decay into e'^e" pair. This method 
is known as internal conversion method [14] and is based on two assumptions. 
The first assumption is that the ratio of direct to inclusive photons is the 
same for real as well as virtiial photons having m~^, < 30 McV, i.c, 7dii./7inci = 
7dir/7inci- Sccondly, the mass distribution follows the KroU-Wada formula [15]: 



(Puee 2a 1 Ami 2m; 



ee 



Sdn^. (2) 



Here, m^ and m^e arc the masses of electron and e^e~ pair respectively and 
a is the fine structure constant. This method is used for Compton scattering 
(g + g — > g + 7* — > g + e"*" + e~), Dalitz decay ( 7r°,r? — > e"'"e~7, w — > e+e~7r°) 
and also for two 7 decay of several other hadrons. The factor S in Eqn. (2) is 
process dependent and for tt" 77* 7e+e~) decay it is expressed as [16]: 

S=\F{ml,)\'^(l-'^J , (3) 

where, Mh is the hadron mass and F{m1^) is the form factor. The factor 
S is for TOee > and goes to 1 as ^ or m^e ^ Pt- The key 
advantage of this method is the greatly improved signal to background ratio 
which is achieved by elimination of the contribution of Dalitz (tt") decay. The 
experimentally measured quantity is the ratio of e+e^ pairs in a particular 
invariant mass bin and the direct photon spectrum is obtained by multiplying 
7dir/7inci ^^ic mcasurcd inclusive photon spectrum (left panel of Fig 3). 

Due to the excellent resolution of the PHENIX detector to measure charged 
particles at low momenta, another powerful technique known as 'tagging of 
decay photons' [17] is used to eliminate the tt" decay background. This method 
is very useful in the low and intermediate pt range ( 1 < < 5 GeV) as 
the systematic uncertainties introduced by detector efficiency, acceptance etc. 
cancel out for measuring j/Jt^o directly in-place of conventional double ratio 
{R = (7/7r°) measured/ (7/7'''') decay) technique. In this method, the invariant 
mass (m^e+e-) distribution of the 7r° ( — > 77* — > 7e+e~) decay products is 
constructed and then mixed event analysis is used for the final subtraction 
(middle and right panel of Fig 3). 
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Fig. 3. Left Panel: Invariant mass distribution of virtual photons from n and 77 
Dalitz decay as well as from direct photons [14]. Middle and right panels: Invariant 
mass spectrum in 'tagging' of decay photon method [17]. 



2.3 Sources of direct photons 



Direct photons can be classified into different categories depending on their 
origin from different stages of the expanding fireball formed after the collision. 
These are: (1) prompt photons, which originate from initial hard scatterings, 

(2) pre-equilibrium photons, produced before the medium gets thermalized, 

(3) thermal photons from quark-gluon plasma as well as by hadronic reactions 
in the hadronic phase, and (4) photons from passage of jets through plasma. 
It is not possible experimentally to distinguish between the different sources. 
Thus, theory can be used with a great advantage to identify these sources of 
direct photons and their relative importance in the spectrum [18]. 



Partonic processes for production of prompt photons 
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(a) (b) 

Fig. 4. Partonic processes for production of photons from (a) quark gluon Compton 
scattering, quark anti-quark annihilation process and (b) quark fragmentation. 



In relativistic heavy ion collisions, prompt photons are produced due to 
quark gluon Compton scattering {q + g g + j), quark anti-quark annihila- 
tion process (9 + 9 — > g + 7), and quark fragmentation ( 9 — > (? + 7) following 
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scattering of partons of the nucleons in the colhding nuclei (see Fig 4). At 
lowest order in aa^, quark gluon Compton scattering and quark anti-quark 
annihilation processes dominate the photon production. In Next to Leading 
Order (NLO) calculation, many more complicated scattering processes ap- 
pear in the photon production cross-section and the total contribution can be 
written as addition of two different terms as [19], 

da da^o) da^P) 

(4) 



dpj-d,ri d p rpdrj d p j^drj 

In the above equation 'D' stands for 'Direct' or the total Compton scattering 
and annihilation contribution while the photons from fragmentation are de- 
noted by 'F'. It is clear that the 'D' photons arc well separated from hadrons. 
The 'F' photons on the other hand have their origin in coUinear fragmenta- 
tion of coloured high px partons and are accompanied by hadrons. This can 
be used with advantage as the produced photons can be separated out by 
employing useful isolation cuts. The two terms in Eqn. (4) can be written 
explicitly as [20]: 

d^^^"^ _ ^ ... .. . 



j dxidx2 Fi/f,^{xuM) Fj/h^{x2,M)- 



^ , ^^K\f\,n^M,M^)\ (5) 



d'^rpdr] 27r 
and 

da(^) ^ r . . dz 



dJ^Tdn ^ _ 
%,3,k=q,q,g 



dxidx2— Fi/hi{xi,M) Fj/h2{x2,M) D^/k{z,Mp) 



Here Fiifi^^{x,M) are the parton distribution functions and as{nR) is the 
strong coupling defined in the MS rcnormalization scheme at the renormal- 
ization scale hr. For details see Ref. [19]. As mentioned earlier, in a complete 
and consistent NLO pQCD [0(Q;as)] calculation important contribution to 
prompt photon result arises from various possible 2^3 process like ab — > ^cd 
for the Direct as well as for the Fragmentation processes [20, 21]. 

In case of photons from fragmentation, the higher order correction is very 
important at low Xt (= '2'PT/^fs)- By choosing a scale for factorization, rcnor- 
malization, and fragmentation, all equal to pt/2, a very good quantitative de- 
scription is obtained for all the available pp and pp data without introduction 
of any intrinsic kr- Results spanning over two orders of magnitude in energy 
and over nine orders of magnitude in cross sections, are shown in Fig. 5. 

It should be mentioned though, that a similar exercise for pions requires 
a scale of pt/3 [22]. 
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In the early calculations, the results for nucleus-nucleus scatterings were 
often obtained by multiplying the pp results for some ^/s with corresponding 
scaling factor iVcou (= (JnnTab] where, unn is the nucleon-nucleon cross- 
section and Tab is the nuclear overlapping function for nuclei A and B) or 
number of binary collisions. In actual practice, often enough, unn was re- 
placed by fjpp. However, as the valence quark structures of protons [uud) and 
neutrons (udd) are different, one needs to correctly account for the iso-spin 
of the nucleons to calculate the prompt contribution. This correction will 
strongly affect the results in the pt range, where the valence quark contribu- 
tion is significant. Also, remember that there is no direct measurement for pn 
and nn cross-sections, though they can be estimated by comparing results of 
scatterings involving deuterons. The effect of shadowing on structure function 
and energy loss of final state quarks before they fragment into hadrons are 
two other corrections [23] which need to be accounted for. 

In order to clearly see this point we show the results for prompt photons 
for nn and pn collisions normalized to those for pp collisions at ^/s = 200 
GeV in the left panel of Fig. 6. It is clear that at low xt (where processes 
involving gluons dominate) the effect is not very significant, however with 
larger values of xt (where the processes involving valence quarks dominate), 
the production of photons decreases by a large factor for the case of pn and 
nn collisions. 
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Fig. 5. World's inclusive and isolated direct photon production cross-sections mea- 
sured in pp and pp collisions compared to JETPHOX NLO predictions, using BFG 
II(CTEQ6M) for fragmentation (structure) functions and a common scale pt/2 [19]. 




Fig. 6. Left panel: Prompt photons for pn and nn collisions normalized to those 
for pp collisions at y/s = 200 GeV. Right panel: Iso-spin, shadowing and energy 
loss corrected Raa for prompt photons at 200A GeV Au + Au collisions using NLO 
pQCD. 
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Fig. 7. Iso-spin, shadowing and energy-loss corrected Raa for prompt photons at 
LHC (left panel) and SPS (right panel) energies using NLO pQCD. 



The transverse momentum dependent nuclear modification factor Raa, 
defined as, 

Raa{pt) = 1 da^^{pT)/dy(PpT 
NcoW da^^{pT)/dy(PpT 

for prompt photons using NLO pQCD and considering iso-spin, shadowing 
and energy-loss effects for 200A GeV Au -\- Au collisions at RHIC is shown 
in the right panel of Fig. 6. We note that, for px < 10 GeV, the iso-spin 
and shadowing corrected result shows an enhancement in the prompt photon 
production compared to the situation when only the iso-spin correction is 
incorporated [24]. This is due to anti-shadowing for larger values of Bjorken 
X or large xt- The inclusion of the energy- loss pushes down the value of Raa 
to less than one for all pr- We also give a comparison with the PHENIX [25] 
experimental data in the bottom lower panel of Fig. 8. 

NLO results at LHC {Pb + Pbm.5A TeV) and SPS {Pb + Pbm7AA 
GeV) are shown in Fig. 7. To clearly demonstrate the relative features of iso- 
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spin, shadowing and energy- loss corrected prompt photons at different collider 
energies, results for Raa a-s function of xt are shown in Fig. 8. 




Fig. 8. Iso-spin, shadowing and fragmentation energy-loss corrected NLO pQCD 
results at different collider energies with different target projectile combination and 
comparison of the RHIC resuhs with PHENIX [25] data. 



2.4 Photon production from Quark-Gluon Plasma 

The thermal emission rate of photons with energy E and momentum p from 
a small system (compared to the photon mean free path) is related to the 
imaginary part of photon self energy by the following relation 

= T^T^M/'^ ^,1 ^ (8) 

where, 77^^''^ is the retarded photon self-energy at a finite temperature T. 
This relation is valid in the perturbative [26] as well as non-perturbative [27] 
limits. It is also valid to all orders in the strong interactions and to order 
in the electromagnetic interactions. If the photon self-energy is approximated 
by carrying out a loop expansion to some finite order, then the formulation of 
Eqn. (8) is equivalent to relativistic kinetic theory. In order to illustrate this, 
we closely follow the treatment of Ref. [28] . 



Electromagnetic probes 11 

Thus using relativistic kinetic theory formulation, the contribution of these 
processes to the rate can be written as [28, 29] : 

where Mi represents the amplitude for one of these processes and the f's 
are the Fermi-Dirac or Bose-Einstein distribution functions as appropriate. 
Positive and negative signs in the last part of Eqn. (9) correspond to Bose 
enhancement and Pauli suppression respectively. 

The integral above can be simplified by introducing Mandelstam variables 
s = (pi + P2)^, t — {pi— PiY and u= {p\— pf'. Now the differential photon 
rate can be written as: 

= (2^1^ / ^"'^^ dE,dE2h{E,)ME,) 

X [1 ± /3(£i +E2- E)]e{Ei +E2- E){aEi^ + bEi + c)-'^/'^{W) 

where, 

b = 2{s + t){Es- E2t), 

c = st{s + t)-{Es + E2tf. (11) 

Considering the photon energy to be large, one can consider, fi{Ei)f2{E2) « 
g-{Ei+E2)/T ^jj^ simpliiy the above as 

^^=(2^?35^'""'V f Ml±.-"-)-/^<|A<.(M)l' (12) 

In the above equation, positive and negative signs stand for fermions (q) and 
bosons {g) in the final state respectively. 

The relation between the amplitude and differential cross-section for mass- 
less particles can be written as, 

da_\Ml_ 

dt ~ Wns^ ^ 

And thus, the differential cross-sections for annihilation process and Compton 
scattering are: 

^^annihilation ^ gTraa, li^ + f 

di ~ ~9s2 ^ ' 

and 
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^^Compton ^ -T^aas V? + 

It ~ ~3^2 ^ ( ^) 

For annihilation process, M = 20 when summing over u and d quarks and 
for Compton scattering Af = 320/3. The total cross-section can be obtained 
after integrating over t. These differential cross sections have a singularity at 
t and/or u = and the total cross-section is infinite as the processes involve 
exchange of mass-less particle. 

To screen this divergence many-body effects are necessary. This approach 
will be discussed later. As a first step let us isolate the region of phase space 
causing the divergences. The integration is done over 

-s + kc^ <t< -kc^, 2fcc^<s<oo, (16) 

where > fc^^ > is an infrared cut-off. 

This treats u and t symmetrically and maintains the identity s + t + u = 
appropriate for all massless particles. 

In the limit that kj^ — > 0, 

l^T'e-^/^[H4ET/k,') + Cp] (17) 
l^T^e-^/^M4ET/k,') + Cb] (18) 
where 

1 12 °° i—W^ 
= TT - CEuier + ^ V In n = 0.0460 , (19) 

n=2 



E 



^■y^Compton 

d^p 

^-^annihilation 



d^p 



Cb = -1- CEuler " ^ V ^ " = -2.1472 (20) 

n=2 

These expressions use the full Fermi-Dirac or Bose-Einstein distribution func- 
tions in the final state. 

These results have a very interesting structure. Thus, the factor 5/9 arises 
from the sum of the squares of the electric charges of the u and d quarks, the 
factor aag comes from the topological structure of the diagrams, a factor 
comes from phase space which gives the overall dimension of the rate, and 
we have the Boltzmann factor e~^^^ for photons of energy E. The logarithm 
arises due to the infrared behavior. 



Infrared contribution 



The infrared divergence in the photon production rate [28] discussed above is 
caused by a diverging differential cross-section when the momentum transfer 
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goes to zero. Often-times long-ranged forces can be screened by many-body 

effects at finite temperatures. Braaten and Pisarski have analyzed problems 
such as this one in QCD [30]. They have argued that a cure can be found in 
reordering perturbation theory by expanding correlation functions in terms of 
effective propagators and vertices instead of bare ones. These effective propa- 
gators and vertices are just the bare ones plus one- loop corrections, with the 
caveat that the one-loop corrections are evaluated in the high temperature 
limit. This makes them relatively simple functions. 

The analysis of Braaten and Pisarski shows that a propagator must be 
dressed if the momentum flowing through it is soft (small compared to T). This 
is because propagation of soft momenta is connected with infrared divergences 
in loops. Dressing of propagators are necessary, otherwise corrections due to 
these are also inflnite. 

Using the one loop corrected propagators and vertices and the contribution 
to the rate coming from the infrared sensitive part of phase space can be 
written as. 




where, 2m,2 = ig^T^. 

Adding this contribution to those given by Eqns. (19) and (20), the final 
result can be written as; 



d^p 9 27r2 \ T J ^ ' 

This is independent of cut-off kc- Thus, the Braaten-Pisarski method has 
worked beautifully to shield the singularity encountered above. We also note 
that in kinetic theory calculation, the effective infrared cutoff is kc^ = 2mq^. 
In earlier works an infrared cut-off was often imposed by giving the exchanged 
quark an effective temperature-dependent mass. 

These early results have been brought to a high degree of sophistication 
and results complete to leading order in with inclusion of LPM effects are 
now available, which should be used for detailed calculation [31]. 



Photons from passage of jets through QGP 

The relativistic heavy ion collisions at RHIC (and LHC) energies are marked 
by a large production of high energy quark and gluon jets which lose energy 
while passing through the QGP due to c;ollision and radiation of gluons. This is 
the celebrated phenomenon of jet- quenching. A quark jet having a transverse 
momentum pt would be formed within r ~ l/pr which can be much smaller 
than To, when QGP is formed, for large pT- This quark (or anti-quark) jet, 
while passing through QGP may annihilate with a thermal anti-quark (or 
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quark) or undergo a Compton scattering with a thermal gluon and lead to 
production of a high energy photon. This process is called the jet-photon 
conversion [9]. 

We have seen that the kinematics of the annihilation of a quark anti- 
quark pair {q + q ^ J + g) is expressed in terms of the Mandclstam variables 
s = {pq + Pq)^, t = {Pg — P-y)^ and u = {j)q — p-yY- We also recall that the 
largest contribution to the production of photons arises from small values of 
t or u, corresponding to p^ ^ pq or p^ ^ pq [2, 9, 32]. 

The phase-space distribution of the quarks and gluons produced in a nu- 
clear collision can be approximately decomposed into two components, a ther- 
mal component /th characterized by a temperature T and a hard component 
/jet given by hard scattering of the partons and limited to transverse momenta 

^ 1 GeV: /(p) = /th(p) + /jot(p)- ./jet dominates for large momenta, while 
at small momenta / is completely given by the thermal part. 

The phase space distribution for the quark jets propagating through the 
QGP is given by the perturbative QCD result for the jet yield [9] : 

^^-(P) = Jq^r^y''^^^'^^ - y^^^^- - - ^) ^''^ 

where gf, = 2 x 3 is the spin and colour degeneracy of the quarks, R± is the 
transverse dimension of the system, and the r] is the space-time rapidity. R{r) 
is a transverse profile function. Tmax is the smaller of the lifetime r/ of the 
QGP and the time taken by the jet produced at position r to reach the 
surface of the plasma. 

One can approximate the invariant photon differential cross sections for 
the annihilation process and Compton scattering as [2, 9], 

~ <7(")(«)^i^7[^(P7 - P,) + '^(P7 - P9-)]> (24) 

and 

^ ~ ^^""^ (s)£^7'^(P7 - Pi) (25) 

Here a'") (s) and a'^^^ (s) are the corresponding total cross sections. 

Using Eqn. (24) and Eqn. (25), the rate of production of photons due to 
annihilation and Compton scattering are given by [2]: 

^^^^ = 2(#g^'(P^) / ^V,(P)[1 + /.(P)].^'^H^) 

x ^^^^'"% (.^^) (26) 
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^'^'-"^'Ki,^,) (27) 



The f's are the distribution functions for the quarks, anti-quarks and gluons. 

Inserting thermal distributions for the ghions and quarks one can obtain an 
analytical expression for these emission rates for an equilibrated medium [2, 
28, 29]. 

In order to estimate the jet photon conversion contribution, wc first note 
that, the integrals over p in Eqns. (26) and (27) are dominated by small 
momenta. Therefore dropping the jet part in the distributions, /(p) in the 
integrands is approximated by the thermal part. Now performing the integrals 
and identifying the quark and anti-quark distributions outside the integrals 
with the jet distributions, results for Compton and annihilation scatterings 
due to jet- conversions are given as. 



'd^p^d^x ^d^p^d^x 

= i^E( V^'t^'^P-) + hi^.)\nH^) + C] (28) 

Here, C = —1.916. If we include three lightest quarkflavours, then^^ e^^^/e^ = 
2/3. We also assume that the mass m introduced here to shield the infrared 
divergence can be identified with the thermal quark mass mth- 

These pioneering works have now been corrected for energy loss and flavour 

change suffered by the jets, as they pass through the plasma [33] as well as 
bremsstrahlung induced by the passage of the jets through the plasma [34]. 



Equilibration time for photons 

It is of interest to get an idea about the equilibration time of photons in the 
medium. Following the treatment of Ref. [28] , we note that the six-dimensional 
phase space distribution dn/d^p{= dN/d^xd^p), for the photons satisfies the 
rate equation: 

d_ / dn \ ^dMf dn/d^p \ 

dt\d^p) d^p\ dn^jd^pj' ^ ' 

In the above equation, dueq/d^p is the equilibrium distribution and is ex- 
pressed as Planck's distribution, 

'^^^'^ -2 1 

W ~ (27r)3 e^/^-l" ^ ' 
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Fig. 9. Typical liadronic reactions for photon production. 



Considering Tcq {— -^/^^) as the time for equilibration and assuming zero 
photons at the beginning we can write, 



d^p d^p 



(l _ e*/^e<.) . (31) 



Using the rate equation, the thermalization time can be expressed in a sim- 
phfied form (considering E > 2T) as: 



9£; 



For energy values E =0.5, 1, 2, 3 GeV, corresponding values for Teq will be 
about 270, 356, 505, 639 fm/c respectively, when the temperature is 200 MeV. 
As the life time of the system is very short, of the order of few tens of fm/c, 
it is very clear that the high energy photons will never reach the equilibrium 
state in heavy ion collisions. This is an important confirmation for the validity 
of the assumption made in all such studies that photons, once produced in 
the collision, leave the system with out any further re-interaction. 



2.5 Photons from hot hadronic matter 

Hot hadronic matter produced after the hadronization of the quark-gluon 
plasma will also lead to production of photons due to hadronic reactions. 
These photons will dominate the spectrum at lower ( < 1 GeV ). The first 
ever calculation of production of thermal photons from hadronic matter was 
performed by Kapusta et al. [28]. 

In a hot hadronic gas (having temperature of the order of pion mass), 
the most important hadronic constituents for photon production are n and 
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p mesons (see Fig. 9). The low mass of pions and the large spin iso-spin 
degeneracy of p mesons, make them the most easily accessible particles in the 
medium. In order to illustrate the photon production from these two mesons, 
we closely follow the treatment of KLS [28]. In a hadronic reaction involving 
charged tt and p meson, the typical Lagrangian describing the interaction can 
be written as: 

C=\D^^ f -ml I # p - Jp^.p"^ + \mlp^ff - ^F^.F"- (33) 

where, 

D^ = d^- leA^ - igpPp., 

# is the complex pion field, p^v is the p field strength and i^^^ is the photon 
field tensor. The diflferential cross-sections for the dominating photon produc- 
ing processes {wp irj) in the hadronic phase is expressed as [28], 




E (GeV) 



Fig. 10. Left Panel: Yield of photons from different hadronic channels (taken 
from [28]). 
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Similarly for 
da 



da , 



dt 



(tt p"*" — > 7r°7) = —{-k'^p 



dt 



s — u — {u — 



7r%) 



".9,; 



48sp: 

2 



2 

cm. 



4(™2 _ 4^2) 



(m — m^Y 



m„ 



t — 



(34) 



^2 



(35) 



and also, 



da Q : , dcr - 

"5^ 



TT 7) 



48sp2 



4(m2 - 4m2)s 
(s - m2 )2 



(s - m2 )2 - Aw?p{rn?p - Ami) 



t 



m„ 



A{ml - Ami) 



s — mz 



s-ml+ ml) 



(36) 



In the above set of equations, s, t, u are the Mandelstam variables and Pc.m. 
is the three momentum of the interacting particles in their centre of mass 
frame. Typical results are shown in Fig. 10. 

Xiong et al. [36] and Song [37] first introduced the 7rp ^ ai — » chan- 
nel for photon production in hadronic phase, whereas baryonic processes and 
medium modification were included by Alam et al. [38, 39]. Several refine- 
ments, e.g., inclusion of strange sector, use of massive Yang-Mills theory and 
t-channel exchange of u mesons, were incorporated by Turbide et al. [35]. The 
last calculation is essentially the state of art result at the moment. Photon 
spectra considering a complete leading rate from QGP [31] and exhaustive 
reactions in hadronic matter [28, 35] are shown in Fig. 11. 
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E (GeV) 

Fig. 11. Complete leading order rates from QGP and exhaustive reactions in 
hadronic matter [28, 31, 35]. 

3 Photons from Pb+Pb@SPS to Au+Au@RHIC 
3.1 SPS 

In order to get an idea of greatly increased insights provided by single pho- 
ton production, let us briefly recall some of the important results from the 
SPS era, many of which preceded the large strides made in our theoretical 
understanding mentioned above. The first hint of single photon production, 
which later turned out to be the upper limit of their production, came from 
the S + Au collisions studied at the SPS energies [40]. 

These results were analyzed in two different scenarios by authors of 
Ref. [41]. In the first scenario, a thermally and chemically equilibrated quark- 
gluon plasma was assumed to be formed at some initial time (tq « 1 fm/c), 
which expanded [42], cooled, and converted into a mixed phase of hadrons 
and QGP at a phase-transition temperature, Tc « 160 MeV. When all the 
quark-matter was converted into a hadronic matter, the hot hadronic gas con- 
tinued to cool and expand, and underwent a freeze-out at a temperature of 
about 140 MeV. The hadronic gas was assumed to consist of tt, p, uj, and 
T] mesons, again in a thermal and chemical equilibrium. This was motivated 
by the fact that the included hadronic reactions involved [28] these mesons. 
This was already a considerable improvement over a gas of mass-less pions 
used in the literature at that time. In the second scenario, the collision was 
assumed to lead to a hot hadronic gas of the same composition. The initial 
temperature was determined by demanding that the entropy of the system 
be determined from the measured particle rapidity density [44]. It was found 
that the scenario which did not involve a formation of QGP led to a much 
larger initial temperature and a production of photons which was consider- 
ably larger than the upper limit of the photon production, and could be ruled 
out. The calculation assuming a quark- hadron phase transition yielded results 
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Fig. 12. Left panel: Upper limits at the 90 % confidence level on the invariant excess 
photon yield for the 7.4 % a^b most central collisions of 200^ GeV ^^S + Au [40]. 
The solid and the dashed curves give the thermal photon production expected from 
hot hadron gas calculation [41], while the dotted curve is the calculated thermal 
photon production expected in the case of a QGP formation. Right panel: Single 
photon production in Pb+Pb collision at CERN SPS. Prompt photons are estimated 
using pQCD (with a K-factor estimated using NLO calculation) and intrinsic kx of 
partons [43]. 

which were consistent with the upper limit of the photon production. These 
results were confirmed [45] by several calculations exploring different models 
of expansion (see left panel of Fig. 12). 

It was soon realized that one may not limit the hadronic gas to contain 
just TT, p, w, and ry mesons, as there was increasing evidence that perhaps all 
the mesons and baryons were being produced in a thermal and chemical equi- 
librium in such collisions. Thus authors of Ref. [46] explored the consequences 
of using a hadronic gas consisting of essentially all the hadrons in the Particle 
Data Book, in a thermal and chemical equilibrium. This led to an interesting 
result for the Pb + Pb collision at SPS energies, for which experiments were 
in progress. It was found that with the rich hadronic gas, the results for the 
production of photons in the phase-transition and no phase transition models 
discussed above were quite similar, suggesting that measurement of photons 
at the SPS energy could perhaps not distinguish between the two cases. How- 
ever, in a very important observation, it was also noted that the calculations 
involving hot hadronic gas at the initial time would lead to hadronic densities 
of several hadrons/fm^, and while those involving a quark gluon plasma in 
the initial state would be free from this malady. Thus, it was concluded that 
the calculations involving a phase transition to QGP offered a more natural 
description. 

The WA98 experiment [11], reported the first observation of direct photons 
in central 158A GeV Pb + Pb collisions studied at the CERN SPS. This was 
explained [47] in terms of formation of quark gluon plasma in the initial state 
(at To ~ 0.2 fm/c), which expanded, cooled and hadronized as in Ref. [46] 
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(see right panel of Fig. 12). An independent confirmation of this approach was 
provided by an accurate description [48] of excess dilepton spectrum measured 
by the NA60 experiment for the same system. 

Once again the results for single photons were analyzed by several authors 
using varying models of expansion as well as rates for production of photons; 
viz., with or with-out medium modification of hadronic properties (see, e.g., 
Ref. [35, 49, 50]). The outcome of all these efforts can be summarized as fol- 
lows: the single photon production in Pb + Pb collisions at SPS energies can 
be described either by assuming a formation of QGP in the initial state or by 
assuming the formation of a hot hadronic gas whose constituents have mas- 
sively modified properties. The later description, however, involved a hadronic 
density of several hadrons/fm'^ , which raises doubts about the applicability of 
a description in terms of hadrons, as suggested by Ref. [46] . 

3.2 RHIC 

The first experimental photon result obtained from the Relativistic Heavy 
Ion Collider was the centrality dependent single photon data for 200^ GeV 
Au + Au collisions by PHENIX collaboration [12], where single photons were 
identified clearly for > 4 GeV. For central collisions, the lower bound on 
photon transverse momenta was even much lower, upto 2 GeV (Fig. 13). This 
data was in good agreement with NLO pQCD results for pp collisions, scaled 
by number of binary collisions ( without considering the iso-spin effect). 

Thermal radiation dominates the direct photon spectrum at lower values 
of pt (< 3 GeV) and the measured slope of the thermal photon spectrum can 
be related to the temperature of the system. 



Fig. 13. Direct photon results for Au + Au collisions (left panel) and p + p collisions 
(right panel) at ^/s^ = 200^ GeV [12]. 
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In Fig. 14 results from several theoretical models based on hydrodynamics 
for thermal photon production at RHIC are shown. The initial temperature 
and thermalization time for Au + Au collision at 200A GeV quoted by several 
theoretical groups [35, 47, 49, 50, 51, 52] are in the range of 450 - 650 MeV 
and about 0.2 fm/c respectively. 

All these calculations are comparable to the experimental data and with 
each others within a factor of 2 and also confirm the dominance of thermal 
radiation in the direct photon spectrum in low and intermediate pr range. 




Fig. 14. Thermal photon spectra for central 200^4 GeV Au-\-Au collision computed 
within different models [52] compared to the expected pQCD prompt photon yield 
and to the experimental total direct photon spectrum measured by PHENIX [25]. 



Indications for jet conversion photons 

We have already discussed the procedure for calculating the high energy pho- 
tons due to passage of jets through the quark gluon plasma. 

The parametrized distribution of jets (quarks, anti-quarks, and gluons) 
obtained by using CTEQ5L parton distribution function and EKS98 nuclear 
modification factor is given by: 

dVd^'^=° - ^■^^d^'^=° - ^^^^ 

Here Taa — 9A^ /8TrR±^ is the nuclear thickness for a head-on collision and 
to include the higher order effects a K factor of 2.5 is introduced. Numerical 
values for the parameters for quarks, anti-quarks and gluons can be found in 
the Ref. [9]. 

Now as a first step, let us ignore the transverse expansion of the plasma and 
assume that a thermally and chemically equilibrated plasma is produced in the 
collision at an initial time tq at temperature Tq. In an isentropic longitudinal 
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expansion, Tg and Tq are related by the observed particle rapidity density 
(dN/dy) by: 



27r4 



1 



dN 



45C(3) 4a7ri?j^ rfy 



(38) 



where a — 42.257r^/90, for QGP consisting of u,d and s quarks and gluons. 
dN/dy can be taken as ~ 1260, based on the charge particle pseudo-rapidity 
density measured by PHOBOS experiment [53] for central collision of Au 
nuclei at Vs^vtv — 200 GeV. For central collision of Pb nuclei at LHC energies 
dN/dy ~ 5625 was used by Fries et al. [9] as in Ref. [54]. 

Further assuming a rapid thermalization, initial conditions can then be 
estimated as [54] Tq = 446 MeV and tq = 0.147 fm/c for the RHIC and 
To = 897 MeV and tq = 0.073 fm/c for the LHC. Now, taking the nuclei 
to be uniform spheres, the transverse profile for initial temperature is given 
by T[r) = To[2(l - rVi?^^)]!/^ where Ri_ = l.2A^/^ fm. The same profile 
R{r) = 2(1 — r'^ /Ri^) is used for the jet production. 

The jet traverses the QGP medium until it reaches the surface or until the 
temperature drops to the transition temperature Tc (~ 160 MeV). 



10 



Photons for Au+Au at s"^=Ax200 GcV 



T from jt'ta through QGP 

Direct photoiis 

— Bremsstrahlung 
Thorraal photons 



10 12 

Pt |GcV| 



14 




7 from jets through QGP 

Direct photons 

— Bremsstrahlung 
Thermal photons 



8 10 12 14 16 18 20 

Pt[GcV: 



Fig. 15. Jet conversion photons at RHIC [left panel] and LHC [right panel] energies 
(taken from Ref. [9]). 



First estimates for the jet conversion photons at RHIC and LHC energies 
along with other sources of photons having large transverse momentum are 
given in Fig. 15. We see that the jet conversion photons make a fairly large 
contribution in the pT range of 4-10 GeV. 

The centrality dependence of jet conversion photons was also studied by 
Fries et al. [55] (see figure 16), which indicates a small but clear contribution of 
photons due to passages of jets through the plasma. These early calculations 
have been brought to a high degree of sophistication by the McGill group [56] 
(see figure 17), where jet quenching and jet-photon conversion along with 
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Fig. 16. Photon yield as a function of pr in central (0-10%) Au + Au collisions at 
200^4 GeV. Results of jet-photon conversion (solid and labeled), primary hard pho- 
tons (dashed) and the sum of the two (upper-most solid curve) along with thermal 
photons are shown separately. Data are from the PHENIX collaboration [25]. 



bremsstrahlung is treated in a single frame work. The details can be seen in 
Ref. [56]. 

In Fig. 17 the results for thermal photons, direct photons due to primary 
processes, bremsstrahlung photons and the photons coming from jets passing 
through the QGP in central collision of gold nuclei at RHIC energies is plotted. 
The quark jets passing through the QGP give rise to a large yield of high- 
energy photons. For RHIC this contribution is dominant source of photons 
up to ~ 6 GeV. Due to multiple scattering suffered by the fragmenting 
partons a suppression of the bremsstrahlung contribution is found. This will 
further enhance the importance of the jet-photon conversion process. 




Fig. 17. Yield of photons in Au + Au collisions at RHIC, for centrality classes 
- 10% (left panel) and - 10% (right panel). See [56] for details. The data sets 
are from Refs. [57] and [25] respectively. 
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Obviously a high statistics data at several centralities and energies will go 
a long way in clearly establishing the presence of jet-conversion photons at 
RHIC and LHC energies. This is also of importance, as there are indications 
that these photons measure the initial spatial anisotropy of the system [56]. 



4 Dileptons 




(a) (b) 

Fig. 18. Dilepton production from (a) quaxk anti-quark annihilation and (b) tt"* 
annihilation through p channel. 



Virtual photons or dileptons are also very powerful and efHcient probes 
like the real photons, to study the dynamics of heavy ion collisions and the 
properties of the medium created in the collisions. Real photons are massless, 
whereas dileptons are massive. Thus, invariant mass M and the transverse 
momentum are the two parameters available for dileptons, which can be 
tuned to investigate the different stages of the expanding fireball. Dileptons 
having large invariant mass and high -pT^ are emitted very early, soon after 
the collision when the temperature of the system is very high. On the other 
hand, those having lower invariant masses come out later from a relatively 
cooler stages and at low temperatures. 

Similar to the real photons, dileptons are also emitted from every stage 
of heavy ion collisions [2]. In the QGP phase, a quark can interact with an 
anti-quark to form a virtual photon 7*, which subsequently decays in to a 
lepton and an anti-lepton pair, together known as dilepton. In the 
hadronic phase, dileptons are produced from interactions of charged hadrons 
with their anti-particles by processes like (7r+ + tt^ ^ /+ + Z^) (see Fig. 18). 
They are also produced from decay of hadronic resonances like p, w, ^, and 
J/tf' as well as from Drell-Yan process. In the Drell-Yan process, a valence 
quark from a nuclcon in the projectile nucleus interacts with a sea anti-quark 
from a nucleon in the target nucleus to form a virtual photon, which decays 
into a lepton pair. The dilepton emissions can be classified into three distinct 
regimes in a rough estimation, depending on the invariant mass M of the 
emitted lepton pairs. These are: 
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Fig. 19. Expected sources of dilepton production as a function of invariant mass in 
relativistic iieavy ion collisions (schematic). 



a) Low Mass Region (LMR): M < M^{= 1.024 MeV): In this mass range, 
vector meson decays are the dominating source of dilepton production and 
medium modified spectral density is one of the key issues which needs to be 
addressed. 

b) Intermediate Mass Region (IMR): < M < Mj/^{= 3.1 GeV): 
In intermediate mass region, continuum radiation from QGP dominates the 
dilepton mass spectrum and thus this region is important for getting a pure 
QGP signature. 

c) High Mass Region (HMR): M > Mj/^: In the HMR, the most interest- 
ing phenomenon are the primordial emission and heavy quarkonia like J/^ 
and T suppression. 

A schematic diagram of continuous dilepton mass spectrum is shown in 
Fig. 19. 

4.1 Dileptons from QGP and hadronic phase 

Dilepton production by quark anti-quark annihilation process or by 7r+ 7r~ 
annihilation process can be expressed in a general form as, 

a++a-^l+ + r (39) 

where, the particle 'a' can be either a quark or a tt meson. We closely follow 
the treatment of Ref. [10] for calculation of dilepton spectra in this section. 
From Quantum Electrodynamics (QED), the cross section for e+e~ — > fJ.^ fJ.'~ 
can be written as, 

47r 

a(M) = ^ 

^ ' 3 M2 



2mf 
M2" 



4m/ 
1P 



(40) 
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In the above equation, a is the fine structure constant (« 1/137), m; is the 
mass of /i, and M is the dilcpton invariant mass. For the quark anti-quark 
annihilation process, this cross-section is modified by inclusion of a colour 
factor Nc{= 3), and the fractional charges of the up(u) and down(d) quarks. 
Thus, the new cross-section takes the form as, 

<7q(M) = F,a(M), (41) 

Fq = 7V,(2s + l)2^e^ (42) 
/ 

where s is the spin of quarks, e/ is the fractional charge, and the sum is over 
different quark flavour 'f. If we consider only u and d quarks having fractional 
charges 1/3 and 2/3 respectively, then the factor Fq is 20/3. 

In the hadronic phase, the tt+tt" annihilation from vector meson domi- 
nance model can be expressed as, 

7r+ -h TT" ^ p ^ Z+ -h r (43) 

For this case, the QED cross-section is multiplied by form factor F„{M) which 
is of Breit-Wigner form as; 

4 

TTl 

Here, Mp is the mass of the p meson (-- 770 MeV) and Fp ( 155 MeV ) is the 
decay width. Thus, the total pion cross-section as a function of M becomes, 

a^{M) = F,a(M)(l - 4ml/M^y/\ (45) 

The reaction rate 'R' can be obtained from the kinetic theory: 

R{a+a~ l~^r) = j ^:^/(Pi) j ^^/(P2)o'(a"^a~ PiP2)wrci 

(46) 

where, 

iilii2 

/(p), is the occupation probability at momentum p and energy E = yj'p^ + rrv^. 
Now, using the distribution function /(p) <^ e~^/'^ , and integrating over five 
out of the six variables, the reaction rate takes a simplified form as, 

R{T) = J2^J^ a{z)z\z' - 4zl)K,{z)dz (48) 

Here, z = M/T, Za = m,a/T, and Ki is the modified Bcsscl function of the 
first kind. The value of the parameter zo is taken as the larger of 2ma/T 
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and 2mi/T. Now, for massless u and d quarks ( as m„, <C T), the e+e 
emission rate takes a simple form of law, 

10 



R- 



97r3 



a 



2rpi 



(49) 



One can also estimate the relaxation time for lepton pairs to come to equilib- 
rium with the QGP: 



Qtt 



For a temperature range of 200 to 500 MeV, the vahie of f^ci varies from 
20 — 60 X 10^ fm/c. We know that the life time of the QGP phase is only 
~ 10 fm/c, which is more than three order of magnitude smaller than tj-ei- 
Thus, like the real photons, the produced lepton pairs also escape the system 
without suffering significant absorption in the medium. 

Now, the dilepton emission rate R is actually defined as the total number 
of lepton pairs emitted from a 4-volume element d'^x (= (Pxt dr] r dr) at a 
particular temperature T is given by: 

R = dN/d^x. (51) 

Thus, the rate of production of dileptons having invariant mass M can be 
expressed using Eqn.(46) as. 



dN 



a{M) 



M^TKi{M/T) 



1 



4m„ 



M2 



(52) 



d^xdM'^ 2(27r)4 

From the last equation and using properties of the modified Besscl's func- 
tions, the production rate per unit 4-v olume for total energy E, momentum 
p and invariant mass M (where, E = y/p'^ + ) can be written as, 

4m^' 



E- 



dN 



a{M) 



-E/T 



1- 



M2 



(53) 



d^xdM'^d^p 4(27r)5' 

Integrating the rate of emission over the entire 4-volume from QGP and 
hadronic phase, one can obtain the pt spectra at a particular M as. 
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fQGp{r,T) 



QGP 



HM 



(54) 



where, rj is the space time rapidity. In the above equation, the temperature 
T, QGP and hadronic matter (HM) distribution functions (/qgpj/hm), all are 
functions of space r (x, y), transverse velocity vt, and proper time r. Eqn. (54) 
can be solved numerically by using any appropriate model [e.g., Ref. [58]] 
with a proper equation of state (EOS). One can also get the invariant mass 
spectrum by integrating out the variable pr from Eqn. (54). The invariant 
mass spectrum of thermal dilepton is dominated by QGP radiation above (j) 
mass and hadronic radiation outshines the QGP contribution for M < Ma,. 
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4.2 Medium modification 

As mentioned earlier, in the low mass region, dilepton emission is largely me- 
diated by /9(770), a broad vector meson, as a result of its strong coupling 
to the TTTT channel and a short lifetime, which is about 1.3 fm/c. In-medium 
properties of vector mesons, like change in medium mass and width have long 
been considered as prime signatures of a hot and dense hadronic medium. 
CERES/NA45 [59, 60] started the pioneering experiment on dilepton mea- 
surement during the period 1989-1992 at CERN SPS. For p + Be and p + Au 
collisions, the theoretical estimates, considering only hadron decay (as source 
of dilepton production), were in excellent agreement with the experimental 
data points. However, the CERES/NA45 experimental data for S + Au and 
Ph + Au at 200^4 GeV and 158 A GeV respectively exhibit excess radiation 
beyond the electromagnetic final state decay of produced hadrons below the 
(j) mass. From these observations, it was concluded that, the theoretical mod- 
els based on tttt annihilation can reproduce the experimental data, only if 
the properties of the intermediate p meson is modified in the medium. This 
was an exciting discovery and a genuine consequence of many body physics. 
However, the resolution and statistical accuracy of the data was insufficient to 
distinguish between models suggesting a drop in the mass of p mesons [61] and 
those which suggest an increase in their decay width [62] and thus, determine 
the in-medium spectral properties of the p meson. 

Shortly after this, an excess dimuon (over the sources expected from pA 
measurement) was identified experimentally by Helios/3 [63] (measured both 
e+e~ and and NA38/NA50 [64] collaborations. 




m„ (GeV/c^) m„ (GeV/c^) 

Fig. 20. Left panel: Invariant mass spectrum of e^e~ pairs omitted in 158^ 
GeV Ph -\- Au collisions from the combined analysis of 1995 and 1996 data by 
CERES/NA45 [59, 60]. The solid line shows expected yield from hadron decay and 
dashed lines indicate the individual contribution to the total yield. Right panel: 
Comparison of the experimental data to i) free hadron decays without p decays 
(thin solid line), ii) model calculations with a vacuum p spectral function (thick 
dashed line), iii) with dropping in- medium p mass (thick dash-dotted line), iv) with 
a medium modified p spectral function (thick solid line) . 
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Rapp and Shuryak [65] and Kvasnikova et al. [48] showed that, the excess 
dimuon observed by NA50 in the mass region 1.5 < M < 2.5 GeV can be 
explained by thermal signal without invoking any anomalous enhancement in 
the charm production. Kvasnikova et al. studied the NA50 intermediate mass 
dimuon result using hydrodynamic model and a detailed analysis of the rates 
of dilepton production from hadronic phase [48]. Detector resolution and ac- 
ceptance were accurately modeled in their calculation and the normalization 
was determined by a fit to the Drell-Yan data using the MRSA parton distri- 
bution function as done in the experimental NA50 analysis. The results are 
shown in Fig. 21. They also studied the centrality dependence of the NA50 
data using hydrodynamics by incorporating azimuthal anisotropy in the cal- 
culation. The results were in fairly good agreement with the measured excess 
dilepton data, as can be seen in the third panel of Fig. 21. 




Fig. 21. (a) The dimuon invariant mass spectrum after correcting for detector 
acceptance and resolution and NA50 data. The Drcll-Yan and thermal contributions 
are shown separately along with correlated charm decay and direct decays of the J/>f' 
and J/^. (b) The dimuon transverse momentum spectrum, (c) Centrality dependent 
results from Ref. [48] and NA50 data [64]. 



Dropping rUp vs. increasing Fp 

The CERES data [59, 60] was unable to distinguish between the models sug- 
gesting drop in the mass of p meson and those suggesting an increase in its 
width, as seen in Fig. 20. The high statistics data obtained by the NA60 Col- 
laboration for In + In collisions at SPS energies for the dimuon excess rather 
clearly ruled out the models advocating the dropping mass scenarios. They 
rather firmly establish the models advocating the substantial broadening of 
the decay width of the p meson due to many body effects (see Fig 22). 

5 Elliptic Flow 

Elliptic flow is one of the key observables in relativistic collisions of heavy nu- 
clei, which confirms the collectivity and early thermalization in the created hot 
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Fig. 22. Left panel: In medium vector meson spectral functions from hadronic many 
body theory. The p in cold nuclear matter at several densities proportional to the 
saturation density po = 0.16 fm"''. Right panel: NA60 [66] excess dimuons in central 
In + In collisions at SPS compared to thermal dimuon radiation using in-medium 
electromagnetic rates [67]. 



and dense matter. For a non-central collision (impact parameter 6 7^ ) of two 
spherical nuclei, the overlapping zone between the nuclei no longer remains cir- 
cular in shape, rather it takes an almond shape. This initial spatial anisotropy 
of the overlapping zone is converted into momentum space anisotropy of par- 
ticle distribution via the action of azimuthally anisotropic pressure gradient, 
which gives rise to elliptic flow [68]. Note that, this anisotropic flow or ellip- 
tic flow can also be produced in central collisions of deformed nuclei, such 
as U + U collisions [69]. The driving force for momentum anisotropy is the 
initial spatial eccentricity (=< — > / < + >) and the momen- 
tum anisotropy can grow as long as e^, > [70] . Elliptic flow coefficient V2 is 
quantified as the second Fourier co-efficient of particle distribution in the pr 
space, which is of the form: 

dN(b) dN(b) , ^ , 

^' - ^' -(6) X [l + 2tii(pT,6)cos(0) 



PTdpT dy d<f) 2npT dpx dy 

+ 2v2{pT, b) cos(2(/)) + 2v3{pT, b) cos(3(/.) + ...]. (55) 

At mid-rapidity [y = 0) and for collisions of identical nuclei, only the even 
cosine terms survive in the above Fourier series and V2 is the lowest non- 
vanishing anisotropic flow co-efficient. Now, the value of V2 depends on the 
impact parameter b, transverse momentum pT as well as on the particle species 
through their rest masses m. For massless real photons, V2 depends on the 
elliptic flow of parent particles. 



5.1 Thermal Photon v-^ 



Elliptic flow co-efficient V2 for photons is a much more powerful tool than the 
V2 of hadrons, to study the evolution history of the system depending on the 
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Fig. 23. Left panel: V2{pt) for thermal photons along with the V2 for tt and p 
mesons. Contributions the from quark matter and the hadronic matter are shown 
separately with the sum contribution. Right panel: Thermal photon V2 for different 
centrality bins [taken from Ref. [7]]. 

different emission time and production mechanism of photons compared to 
hadrons. Photons are emitted from all stages and throughout the evolution of 
the system, whereas, hadrons are emitted only from the freeze-out surface at a 
relatively much cooler temperature (~ 100 MeV). The interplay of the photon 
contributions from fluid elements at different temperatures with varying radial 
flow pattern, gives the photon V2 a richness, which is not possible for hadrons. 
Also, photons emitted from QGP phase as a result of qq annihilation and 
quark gluon Compton scattering, carry the momenta of the parent quarks or 
anti-quarks, which makes them an unique tool for QGP study. 

Ideal hydrodynamic model has successfully predicted the elliptic flow 
for different hadronic species upto a pt value of 1.5(2.3) GeV for mesons 
(baryons) [70] at RHIC. Experimental data show that, V2{pt) for hadrons sat- 
urates beyond that px range. However, ideal hydrodynamics predicted 1^2 (pt) 
rises monotonically with px- The saturation at higher px is is often explained 
as due to viscosity, neglected in the discussion here. 

Fig. 23 shows the pr dependent V2 for thermal photons [7] for a typical 
impact parameter b = 7 fm. Contributions from quark matter (QM) and 
hadronic matter (HM) along with the sum of the two are shown separately 
(left panel of Fig. 23). One can see that, V2{HM) rises monotonically with 
Pt, similar to the hadronic V2 predicted by hydrodynamics. On the other 
hand, V2{QM) is very small at high px or early times, as very little flow 
is generated by that time. V2{QM) rises for smaller values oi pt and after 
attaining a peak value around 1.5 — 2.0 GeV, it tends to as — > 0. The 
total flow or W2(QM-|-HM) tracks the V2{QM) at high pT in-spite of very 
large V2{HM), as the yield of photons from hadronic phase is very small at 
high Pt- It is well known and mentioned earlier that, for photon energy larger 
than the rest masses of the photon emitting particles, the photon production 
cross-section peaks when the photon momentum and momentum of photon 
emitting particles become almost identical [2]. Thus, at high pT photon V2 





Fig. 24. Left panel: pr integrated elliptic flow of thermal photons from different 
phases at RHIC. Right panel: Time evolution of elliptic flow from different phases, 
normalized by the final value at the time of freeze-out [taken from Ref. [71]]. 



reflects the anisotropies of the quarks and anti-quarks at early times. Also, 
the collision induced conversion of vector mesons (p) starts dominating the 
photon production for px > 0.4 GeV, which gives rise to a structure at the 
transition point in the photon V2 curve [7]. As HM contribution dominates 
the pt spectrum and hydrodynamics is well applicable around that pt range, 
the structure is expected to survive in the experimental result also. 

Recently Turbide et al. [56] have shown that, total contribution to pho- 
ton V2 from prompt fragmentation (small -I- ve V2) and jets-conversion (small 
— ve V2) photons is very small, almost equal to zero. Also, we know that the 
prompt photons from Compton and annihilation processes do not contribute 
to elliptic flow as their emission is not subjected to collectivity and has az- 
imuthal symmetry. Thus, at low and intermediate pT range, V2 from thermal 
photons plays a dominant role in deciding the nature of the direct photon V2- 

For central collision of two spherical nuclei, the spatial eccentricity (e) of 
the overlapping zone is zero, thus the flow co-efflcient V2 is also zero. The 
V2{pt) for thermal photons for different centrality bins are shown in the right 
panel of Fig. 23, where the value of V2 rises from central towards peripheral 
collisions. Now, the elliptic flow as well as spatial eccentricity both rise with 
higher values of impact parameter. However, the ratio of the two remains 
independent of impact parameter up to a very large value of b [71]. The 
Pt integrated V2 for thermal photons from different phases as a function of 
collision centrality shows an interesting behavior (left panel of Fig. 24). The 
V2{h) rises with b until for very peripheral coUisions. For these the system size 
itself becomes very small to generate enough pressure gradient and elliptic 
flow, and as a result V2 {h) decreases. The pr dependent as well as pr integrated 
time evolution results [71] of thermal photon spectra and elliptic flow show 
explicitly the gradual build up of spectra and V2 with time very well. The 
photon V2 from QGP phase saturates within about 4-5 fm/c (right panel of 
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Fig. 24), whereas V2 from hadronic phase is not very significant at early times 
and it saturates much later. 

5.2 Thermal Dilepton V2 

Elliptic flow of thermal dileptons is another very interesting and illustrative 
observable which gives information of the different stages of heavy ion colli- 
sion depending on invariant mass M and transverse momentum pT [72]. At 
the p and (f) masses, the pr spectra and V2{pt) show complete dominance 
of hadronic phase. At these resonance masses radiation from quark matter 
becomes significant only for very large pt{> 4) GeV (left panel of Fig. 25). 
The V2{pt) &t M — nip also shows similar nature as spectra and the total 
V2 is almost similar to hadronic V2 upto a large px- Right panel of Fig. 25 
shows dilepton V2 at rho mass along with V2 of rho meson where, V2{HM) 
tracks V2{p). Also, the hadronic V2 is little smaller than the V2{HM) for dilep- 
tons. The effective temperature of dilepton emission is a little larger than for 
hadrons, thus the spectra of the later is boosted by a somewhat larger radial 
flow [72], which results in smaller elliptic flow. 




Pt (GeV) Pt (GeV) 

Fig. 25. Pt spectra [let panel] and V2 [right panel] for thermal dileptons at M = nip. 
V2{p) is also plotted in the V2 curve for comparison [72]. 

A totally different situation emerges for the pt spectrum and V2 for dilep- 
ton invariant mass M > M^. For M = 2 GeV, the dilepton V2 is similar to 
V2{QM) inspite of the 20 times larger V2{HM), as hadronic dileptons are very 
few compared to QGP radiation at this M value (Fig. 26). 

The Pt integrated spectra and V2 as a function of invariant mass M show 
well defined peaks at the resonance masses (p, lo and 0), as can be seen 
Fig. 27. For M < M^f, the dilepton spectrum is totally dominated by hadronic 
phase and above that it is dominated by contribution from QGP. V2{QM) is 
very small at large M and rises for smaller M. It shows a nature similar to 
the thermal photon V2- The M dependent V2{HM) is significant only at the 
resonance masses and beyond (j) mass its contribution to total V2 is negligible. 
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26. Pt spectra [let panel] and 112 [right panel] for thermal dileptons at M = 2 
[72]. 



Thus, for large values of M or at early times flow comes from QGP phase 
and its value is very small, whereas at later times or at low M values it is 
from HM, which is very strong. The measurement of flow parameter at high 
M and/or high px values can be very useful to reveal a pure QGP signature. 




M (GeV) M (GeV) 

Fig. 27. Left panel: The mass spectra of thermal dileptons from a hydrodynamical 
simulation of central 200^ GeV Au + Au collision. Quark and hadronic matter con- 
tributions are shown separately. Right panel: pr integrated elliptic flow parameter 
for dileptons and various hadrons [72]. 



5.3 Elliptic flow of decay photons 

We have already discussed that most of the produced photons in heavy ion 
collisions are from the 2-7 decay of tt*^ and 77 mesons and the direct photons 
contain a small fraction of the inclusive photon spectrum. Thus it is very 
interesting to know the nature of V2 of decay photons from pion and rj decays. 
The momentum distribution of decay photons from tt^ decay in an invariant 
form can be expressed as [73], 
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fco^(p, fc) = ^'5(p-fc- im^), (56) 
where p and k are the 4-momentum of the pion and photons and m is the pion 
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Fig. 28. Left panel: Elliptic flow of photons from tt" decay at mid-rapidity, /(pt) 
stands for the momentum distribution of the tt". Right panel: Elliptic flow param- 
eters for photons from decay of vr" obtained using the recombination model [taken 
from Ref. [74]]. 



mass. Thus the Lorentz invariant cross-section of the decay photons using the 
decay kinematics from the Eqn. (56) is, 

da f (Pp / da\ 1 1 

Layek et al. [74] have calculated elliptic flow of decay photons considering sev- 
eral azimuthally asymmetric pion distributions (Fig. 28). They have shown 
that fc-r dependent V2 of decay photons closely follows the V2{pt) of n° evalu- 
ated a,t Pt ^ k-T + S (where 6 ^ 0.1 — 0.2 GeV). Similar results were obtained 
for decay photons from r] mesons also. This study could be useful in identify- 
ing additional sources of photons as the V2 of 7r° is similar to that of 7r+ and 
7r~. Also by using the property of quark number scaling or the recombination 
model, the decay photon V2 can help to estimate the V2 of constituent partons 
in the tt" or rj mesons (right panel of Fig. 28). See Ref. [74] for details. 

Experimental measurement of direct photon V2 

PHENIX has measured direct photons and its V2 by subtracting V2 of decay 
photons (2-7 decay of tt" and 77 mesons) from inclusive photon V2 using appro- 
priate weight factor [75]. The procedure followed by them to estimate direct 
photon V2 is as follows: 
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where R is the direct photon excess over hadron decays defined as, 

(7/7r°)bkgd. 

The factor R is measured from spectral analysis. The inclusive and hadron 
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Fig. 29. Elliptic Flow of direct photons from 200^4 GeV Au + Au collision at 
RHIC [76]. 



decay V2 are measured using reaction plane method where the background 
photons are measured using a Monte Carlo simulation. We have already seen 
that the photon V2 shows different signs and/or magnitude depending on the 
production procedure of photons in heavy ion collisions. Thus experimental 
measurement of photon V2 can be a very powerful tool to disentangle the 
intermix contributions from different sources in different px ranges. However, 
the experimentally measured preliminary PHENfX data still contain a large 
systematic error (Fig. 29), which need to be reduced before making any specific 
conclusion about direct photon V2. 



6 Photon tagged jets 



Jets are hard phenomenon characterized by a large momentum transfer be- 
tween partons and are characterized by several hadrons in a small angle around 
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a leading particle as a result of jet fragmentation. The hadrons surrounding 

a jet are known as associated particles of the jet. 

One of the key observables at RHIC energy is the strong suppression of 
leading hadron yield at high values of transverse momentum compared to 
expectations based on p + p or d + A collisions at same collision energy. This 
is the celebrated phenomenon of 'jet quenching'. It is assumed that the hard 
partons lose a large fraction of their energy when passing through the strongly 
interacting medium. A variety of qualitatively different models are available 
using collisional energy loss and radiative energy loss mechanisms to describe 
the nuclear modification factor Raa, defined as, 



Here, Taa is the nuclear overlap function and aNN is the nucleon-nucleon 
cross-section. The major problem for studying jet properties arises from the 
fact that, conventional calorimetric study can not measure the jet energy loss 
very accurately. Also it is very difficult to directly measure the modification 
of jet fragmentation function and jet production cross-section. 

The 'jet quenching' in heavy ion collisions can be studied by measuring the 
Pt distributions of charged hadrons in opposite direction of a tagged direct 
photon. Quark-gluon Compton scattering and quark anti-quark annihilation 
process are the main mechanisms by which direct photons are produced at very 
high Pt, and jets are produced in the opposite direction of these photons. Prom 
momentum conservation law, the initial transverse energy of the produced jets 
are equal to that of the produced photons, i.e, = £jet- Thus uncertainties 
regarding the jet cross-section can be avoided by tagging a direct photon in 
the direction opposite to the jet (see. Ref. [77]). 

Medium modified parton fragmentation function Di^i^^z) ( z is fractional 
momenta of the hadrons) is used to study the jet energy loss. The differential 
Pt distribution of hadrons from jet fragmentation in the kinematical region 
{A(l),Ay) using D^/aiz) can be written as: 



where, C{AyA(p) = j\y\<Ay/2dy !\4>-^^\<A^/2d(t>f{.y -\(t>- is an overall 



factor and }'{}), (t>) is the hadron profile aroimd the jet axis and ra{E^) is 
the fractional production cross section of a typical jet associated with the 
direct photons. Thus, comparison between extracted fragmentation function 
in AA and pp collisions can be used to determine the jet energy loss as well as 
the interaction mean free path in the dense matter produced in high energy 
heavy ion collisions. Thorston Renk [78] has shown that by gamma hadron 
correlation measurement, the averaged probability distributions for quarks 
are accessible experimentally and he has also explained an analysis procedure 



RAA{pT,y) = 



(fNAA/dpTdy 



(60) 



TAAm-'a^^/dpTdy 




(61) 
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capable of distinguishing between different energy loss scenarios leading to 
the same nuclear suppression factor. 
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Fig. 30. Left panel: Nuclear suppression factor Raa for different toy models and 
hydrodynamical simulation compared with PHENIX data. Right panel: Momentum 
spectrum of hard hadrons correlated back to back with a photon trigger normalized 
to the expectation of geometrical absorption (taken from Ref [78]). 



Isolating the bremsstrahlung photons 

In order to measure direct photon cross section by suppressing the hadron 
decay background or in particular to isolate bremsstrahlung photons from ac- 
companying hadrons, an interesting method known as 'isolation cut' can be 
used successfully. The basic assumption for performing the cut is that the 
hadronic energy in a cone around the photon is less than a certain fraction of 
the photon energy, i.e, a photon is considered as isolated if the combined en- 
ergy of the accompanying hadrons is less than eE^ ( where, is the photon 
energy) inside a cone of half opening angle 6 around the photon. The parame- 
ter e is very small ( ^ 0.1) and is called the energy resolution parameter. The 
cone around the tagged photon is known as isolation cone. The cone opening 
can be related to the radius i? of a circle centered around the photon in the 
centre of mass system, where R is defined in terms of rapidity rj and azimuthal 
angle </> as, 

R > \J Ajf + A(j? (62) 

For small rapidities, the half opening angle 5 equals the radius R. This method 
was successfully implemented into the theoretical study of isolated prompt 
photon production considering fragmentation contribution also in next-to- 
leading order (NLO) by Gordon and Vogelsang [21]. A very good accuracy 
of this method over a wide range of isolation parameters for prompt photon 
production was demonstrated in their calculation. Results from RHIC and 
LHC for this would be very valuable. 
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Fig. 31. Left panel: The production rate (per event) of hard photons in a central 
collision of gold nuclei at ^/sNN = 200 GeV as a function of time in the centre of 
mass system. Right panel: Spectrum of photons from various sources. QM and HM 
denote quark matter and hadronic matter contribution respectively (figures are from 
Ref. [6]. 

Dilepton tagged jets 

In the study of photon tagged jets, the main problem arises from the jet 
pair production background, where a leading tt*^ in the jet is misidentified 
as a photon. For an event having a huge background 'isolation cut' is not 
a very efficient mechanism in the low px range to study photon tagged jets 
in heavy ion collisions and it is useful only for pp collisions. In the high px 
range although the background contribution related problems are reduced, 
a substantial problem in photon isolation is created by small opening angle. 
On the contrary jets tagged by dileptons are not affected by background and 
can be used to observe px imbalance, a signal of medium induced partonic 
energy loss. As we have mentioned earlier, for dileptons not only px but also 
invariant mass M is another equally important parameter which can be used 
accordingly to study dilepton tagged jets in the medium. At high px and high 
M, the dilepton yield is much lower compared to low px and low M range, 
however the relative background contribution is also lower in that range. 

At very high transverse momentum Drell-Yan process (ft.i+/i2 ^ l^l^ +X) 
dominates the dilepton production from QGP phase. Srivastava et al. [79] have 
estimated the results for dilepton tagged jets by studying Drell-Yan process 
at NLO in relativistic heavy ion collisions at RHIC and LHC energies. They 
have also shown that correlated charm and bottom decay are unimportant as 
background for dileptons having large transverse momentum or in the kine- 
matical region of interest for jet quenching. Lokhtin et al. [80] have studied 
the dimuon + jet production (including both j* /Z — > ^~ modes) at LHC 
energy and have shown the px imbalance between jX^ n~ pair and a leading 
particle in a jet is clearly visible even for moderate energy loss. It is directly 
related to absolute value of partonic energy loss and almost insensitive to the 
angular spectrum of emitted gluons and to experimental jet energy resolution. 
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Fig. 32. Left panel: The outward, side-ward and longitudinal correlations of direct 
photon predicted by PCM at Kt = 2 GeV, inclusion of thermal photons changes 
the results only marginally. Right panel: Intensity correlation of photons at 1 GeV, 
considering only PCM(BMS), only thermal (hyd), and all PCM-|- thermal photons 
(BMS-Hhyd) (see Ref. [6]). 



7 Intensity interferometry of thermal photons 

The quantum statistical interference between identical particles emitted in 
relativistic heavy ion collisions, provide valuable insight about the shape and 
size of the particle emitting source. We know that direct photons emitted 
from different stages of the collision dominate the pr spectra depending on 
the range of transverse momentum range. Thus, one can extract space time 
dimension of the system at different stages of the collision by measuring the 
correlation radii for photons at different transverse momentum. 

Two particle correlation function C{q, K) for photons having momenta 
fci and ^2, emitted from a completely chaotic source can be written as, 

In the above equation the factor 1/2 appears for averaging over spin and 
S{x, K) is the space time density function defined as. 



(64) 
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where, 

q = k^-k2,K = {k^+k^)/2 (65) 

The correlation function C{q, K) can be expressed in terms of outward, side- 
ward and longitudinal momentum difference and radii ( gout) Q'side; ?iong and 
-Rout, -Rside, -Riong respectively )as, 

C(9out, 9side, 9long) ~ 1 + ^ exp[-(g^^t-Rout + ^side-^side + 9lo„g^Lig)/2] (66) 

Photon four momentum in terms of transverse momentum Ut, rapidity y and 
azimuthal angle can be expressed as 

= (fcr cosh y, fcT cos ■(/;, fcr sin f/;, /ct sinhy) (67) 

and 

giong = \kiz- k2z \ = |fciTsinh?/i - A;2Tsinh?/2| 
Q'out = qr-^T/kr, 
gsido = Iqt - goutkx/^Tl 

Photon interferometry in the QGP phase has been investigated by several 
theoretical groups and the first experimental results on direct photon was 
obtained by WA98 [81] collaboration. Bass, Miiller, Srivastava [6], have cal- 
culated two body quantum correlation of high energy photon using Parton 
Cascade Model (PCM) and ideal hydrodynamic model for central 200A GeV 
Au + Au collision at RHIC. They have shown that one can differentiate be- 
tween the direct photons from early per-equilibrium stage and the same from 
later QGP and hadronic gas stages depending on features of the correlation 
function. Left panel of Fig. 31 shows that about 88% of the total photons hav- 
ing transverse momenta (1 ± 0.1) GeV (produced by hard parton scattering 
in PCM [82]) are emitted within a time period of 0.3 fm/c at rapidity y = 0. 
Right panel of the same figure shows that at fcr = 1 GeV, the contribution 
of pre-equilibrium phase and thermal photons are similar to the total photon 
yield, whereas at fcr = 2 GeV hard photons from PCM outshine the thermal 
ones by an order of magnitude. The emission time in PCM is very small and 
hydrodynamic calculation with tq = 0.3 fm/c allows a smooth continuation 
of the emission rate. Intensity correlation for A;t > 2 GeV (Fig. 31) reveals 
a pre-thcrmal photon dominated small size source of brief duration. On the 
contrary, for fc^ = 1 GeV we can see much larger radii for an extended source 
and suppression of pre-thermal contribution over thermal. 

Intensity interferometry of thermal photons having transverse momentum 
kr ~ 0.1 — 2.0 GeV [43] provide an accurate information about the temporal 
and spatial structure of the interacting medium. 

In reference [43] WA98 data is compared with theoretical results and pre- 
diction are given for RHIC and LHC energies. One dimensional correlation 
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Fig. 33. One dimensional correlation function for the kinematic window used in 
WA98 experiment, assuming a fully source and emitting only single photon. Out- 
ward side-ward and longitudinal correlation function for thermal photons in central 
collision at RHIC energy, (taken from Ref. [43]). 

function in terms of Qinv (invariant momentum difference) for different kx and 
rapidity windows along with WA98 data is shown in figure 33. Theoretical re- 
sults are well fitted in the form C = f + 0.5 aexp[—qf^^Rf^^/2] and are in 
reasonable agreement with the experimental data. At RHIC energy, contribu- 
tion from quark matter increases and as a result the two source aspect in the 
outward correlation radii become more clear. Similar results are obtained for 
LHC energies. It is found that the transverse momentum dependence of the 
different radii are quite different from the corresponding results for pions and 
do not decrease as l/y^niT- For SPS, RHIC, and LHC energies the longitu- 
dinal correlations show similar values, which can be explained as a result of 
boost invariance of the flow pattern. 

8 Epilogue 

We have tried to give a reasonably complete introduction to the exciting 
possibilities provided by radiation of direct photons and dileptons for the 
study of the dynamics of relativistic heavy ion collisions. 
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While low mass dileptons provide insights into the medium modification of 
vector mesons, those having intermediate masses carry signatures of thermal 
radiation from the quark gluon plasma. We have, due to lack of space, left 
out the discussion of correlated decay of charm and bottom mesons which 
give a large contribution to dileptons. These arc; also important as they carry 
valuable information about the extent of thermalization, elliptic flow, and 
energy- loss exhibited by heavy quarks. 

We have discussed sources of direct photons, consisting of prompt photons, 
thermal radiation from quark and hadronic matter and those due to passage 
of jets through QGP. We have also seen that photons carry information about 
the initial temperature, evolution of elliptic flow, and size of the system. High 
statistics data at RHIC and LHC along with photon or dilepton tagged jets 
will go a long way in seeing that we realize the full potential of electromagnetic 
probes of quark gluon plasma. 
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